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Abstract 

We estimate the leading order amplitudes for exclusive photon and me- 
son electroproduction reactions at large Q 2 in the valence region in terms 
of skewed quark distributions. As experimental investigations can currently 
only be envisaged at moderate values of Q 2 , we estimate power corrections 
due to the intrinsic transverse momentum of the partons in the meson wave- 
function and in the nucleon. To this aim the skewed parton distribution 
formalism is generalized so as to include the parton intrinsic transverse mo- 
mentum dependence. Furthermore, for the meson electroproduction reactions, 
we calculate the soft overlap type contributions and compare with the leading 
order amplitudes. We give first estimates for these different power corrections 
in kinematics which are relevant for experiments in the near future. 
PACS : 12.38. Bx, 13.60.Le, 13.60.Fz, 13.60.Hb 

I. INTRODUCTION 

Much of the internal structure of the nucleon has been revealed during the last two 
decades through the inclusive scattering of high energy leptons on the nucleon in the Bjorken 
-or "Deep Inelastic Scattering" (DIS)- regime defined by Q 2 \v — > oo and xb = Q 2 jlMv 
finite. Simple theoretical interpretations of the experimental results can be reached in the 
framework of QCD, when one sums over all the possible hadronic final states. 

With the advent of the new generation of high-energy, high-luminosity lepton accelera- 
tors combined with large acceptance spectrometers, a wide variety of exclusive processes in 
the Bjorken regime are considered as experimentally accessible. In recent years, a unified 
theoretical description of such processes has emerged through the formalism introducing new 
generalized parton distributions, the so-called 'Off- Forward Parton Distributions' (OFPD's), 
commonly also denoted as skewed parton distributions. It has been shown that these dis- 
tributions, which parametrize the structure of the nucleon, allow to describe, in leading 
order perturbative QCD (PQCD), various exclusive processes in the near forward direction. 
The most promising ones are deeply Virtual Compton Scattering (DVCS) and longitudinal 
electroproduction of vector or pseudoscalar mesons (see Refs. [0- || and references therein 
to other existing literature). Maybe most of the recent interest and activity in this field has 
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been triggered by the observation of Ji [|IJ, who showed that the second moment of these 
OFPD's measures the contribution of the spin and orbital momentum of the quarks to the 
nucleon spin. Clearly this may shed a new light on the "spin-puzzle" . 

In section |J we introduce the definitions and conventions of the skewed parton distri- 
butions. We also present the modelizations for these distributions that will be used in our 
cross section estimates. 

In section [TTT| , the leading order PQCD amplitudes for DVCS and longitudinal electro- 
production of mesons are presented in some detail. This is an extension of our previous 
works 

In section [TV], we investigate the power corrections to the leading order amplitudes when 
the virtuality of the photon is in the region Q 2 m 1 — 20 GeV 2 . The correction due to 
the intrinsic transverse momentum dependence is known to be important to get a successful 
description of the 7r°7*7 transition form factor, for which data exist in the range Q 2 ~ 1 — 10 
GeV 2 . For the pion electromagnetic form factor in the transition region before asympto- 
tia is reached, the power corrections due to both the transverse momentum dependence 
and the soft overlap mechanism are quantitatively important. We will therefore take these 
form factors as a guide to calculate the corrections due to the parton's intrinsic transverse 
momentum dependence in the DVCS and hard meson electroproduction amplitudes. In ad- 
dition, for the meson electroproduction amplitude, an estimate is given for the competing 
soft overlap mechanism, which - in contrast to the leading order perturbative mechanism - 
does not proceed through one-gluon exchange. 

In section [V[ we give numerical estimates for different observables for DVCS and for 
the electroproduction of vector mesons p < - ±,0 ' ) , uj and pseudoscalar mesons n^' ^, r\. We give 
several examples of experimental opportunities to access the OFPD's at the current high- 
energy lepton facilities : JLab (E e > 6 GeV) J7|,§, HERMES (E e =27 GeV) and COMPASS 
(£ M =200 GeV). 

Finally, we give our conclusions in section [VI[ 

II. REVIEW OF SKEWED PARTON DISTRIBUTIONS 

A. Preliminary : bilocal operators and gauge invariance 

In the following we often encounter bilocal products of quark fields ip at equal light-cone 
time of the form 

4{-y/2)if>(y/2) \ y+=0 , (1) 

where y + is the light-cone time, with y^ 1 = (y° ± y 3 ) / y/2. At leading twist one has the 
further restriction y± = (y 1 ^ 2 , 0) = 0, that is the bilocal product needs only be evaluated 
along a light-cone segment of length y~. For our purposes we need to consider the more 
general case where only y + is set to zero ( y + = is understood in the remaining of this 
section). 

The product of fields at different points is not invariant under a local color gauge trans- 
formation. To enforce gauge invariance, one replaces Eq. (P by 

^(-y/2)L(-y/2 - y/2)^(y/2), (2) 
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where the link operator L is defined by 

L{-y/2^y/2) = P{e- i So^). (3) 

In Eq. (^|), the path ordered product is evaluated along a curve C joining the points — y/2 
and y/2 and A is the matrix valued color gauge field. So in a general gauge, bilocal products 
of quark fields are not defined independently of the gauge field. To simplify the analysis it 
is therefore convenient to choose a curve C and a gauge such that the link operator reduces 
to the identity. 

Among the many possibilities we choose the curve shown on Fig. |], where it is understood 
that y + = 0. We then split L into 3 factors Li,L 2 ,L 3 corresponding to the segments 
Ci,C 2 ,C 3 indicated on Fig. |[ We first go to the gauge where A + = everywhere. This 
is always possible if one ignores the problems associated with the boundary conditions at 
infinity. Since we only deal with localised systems this is not a true restriction. In this 
gauge we have L\ = L% — 1. While staying in the gauge A + = we still have the freedom 
to make a gauge transformation which depends only on (y 1 ,?/ 2 ). This allows us to go to a 
gauge where y l A l + y 2 A 2 = in the plane y~ = 0. This is the 2-dimensional version of the 
radial, or Schwinger-Fock, gauge ||. In this gauge we have L 2 = 1 which completes our 
argument, that is there exists a curve and a gauge such that the link operator is the identity 
operator. In the following it will always be understood that the bilocal quark field products 
(at y + = 0) are in this gauge [], which is thus a convenient choice for model evaluations of 
the matrix elements of these operators. 



B. Definitions of skewed parton distributions 

To set the framework of this work, we briefly review how the nonperturbative nucleon 
structure information enters the leading order PQCD amplitude for DVCS and hard meson 
electroproduction. 

For DVCS, Ji and Radyushkin f2| have shown that the leading order amplitude in the 
forward direction can be factorized in a hard scattering part (which is exactly calculable in 
PQCD) and a soft, nonperturbative nucleon structure part as is illustrated in Fig. §. In these 
so-called "handbag" diagrams, the nucleon structure information can be parametrized, at 
leading order PQCD, in terms of four generalized structure functions, which conserve quark 
helicity. It was shown in Refs. []ID,ffT , that at leading twist, two more functions appear which 



involve a quark helicity flip. However, the authors of Ref. |T2| have shown that in reactions 
which would involve these helicity flip distributions (such as exclusive electroproduction of 
transversely polarized vector mesons), this contribution vanishes due to angular momentum 



and chirality conservation in the hard scattering. This was argued |L2| to hold at leading 
order in 1/Q and to all orders in perturbation theory. In this work, we will restrict ourselves 
to the four distributions which conserve quark helicity. Then, the matrix element of the 
bilocal quark operator, representing the lower blob in Figs. 0a and 0b, can be expressed at 
leading twist, in the notation of Ji, in terms of the OFPD's H, H, E, E as : 



or any other one which allows to replace the link operator by the identity 
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where ip is the quark field, N the nucleoli spinor and m^ the nucleon mass. In writing down 
Eq. (H) one uses a frame where the virtual photon momentum and the average nucleon 
momentum P M (see Fig. ^| for the kinematics) are collinear along the z-axis and in opposite 
direction. We denote the lightlike vectors along the positive and negative z-directions as 
p» = P+/V2(l, 0, 0, 1) and n" = 1/P+ ■ 1/^(1, 0, 0, -1) respectively, and define light-cone 
components a ± by a 1 * 1 = (a ± a 3 )/y / 2- In this frame, the physical momenta have the 
following decomposition : 



5 (P M + P**) 
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(5) 

(6) 
(7) 
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where the variables m , £ and £ are given by 



m 
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(9) 
(10) 
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In Eq. (H), the OFPD's if 9 , P 9 , H q , E q are defined for one quark flavor (q = u,d and s). The 
functions H and P are helicity averaged whereas H and P are helicity dependent functions. 
They depend upon three variables : x, £, t. The light-cone momentum fraction x is defined 
by k + = xP + . The variable £ is defined by A + = — 2£ P + , where A = p — p and i = A 2 . 
Note that 2£ — > xb/(1 — xb/2) in the Bjorken limit. The support in x of the OFPD's is 
[—1,1] and a negative momentum fraction corresponds with the antiquark contribution. 

To simplify the presentation, we will take £ = £ in the following. As shown in Eqs. (||- 
p"T|) , these two variables have the same value in the Bjorken limit. The kinematic variable 
£, which represents the longitudinal momentum fraction of the transfer A, is bounded by 



n V /Z A72 
< £ < — < 1 . 



m 



(12) 
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A glance at Figs. |2| a,b shows that the active quark with momentum k — A/2 has longi- 
tudinal (+ component) momentum fraction x + whereas the one with momentum k + A/2 
has longitudinal momentum fraction x — £. As noted by Radyushkin ||, since negative 
momentum fractions correspond to antiquarks, one can identify two regions according to 
whether \x\ > £ or \x\ < £ . 

• When x > £, both quark propagators represent quarks, whereas for x < — £ both 
represent antiquarks. In these regions, the OFPD's are the generalizations of the 
usual parton distributions from DIS. Actually, in the forward direction, the OFPD's 
H and H respectively reduce to the quark density distribution q(x) and the quark 
helicity distribution Aq(x) : 

H q (x,0,0) = q(x) , H q (x,0,0) = Aq(x) , (13) 



where q and Aq are defined by the Fourier integrals on the light cone [113 



q ( x ) = V — I (hr< -r « <,,|r(0)-,./»'(/7)|/'> 



v + =y±=o 



Aq(x) = £ / dy~ eixp+v ~ (pS\\ l'0(O)'7- n T5'0(z/) \pS\\) 



(14) 
(15) 



y + =y 1=0 



In Eqs. (0, |15D, p represents the initial nucleon momentum in the DIS process and S\\ 
is the longitudinal nucleon spin projection. From Eqs. Ql4] , |T3J) the quark distributions 
for negative momentum fractions are related to the antiquark distributions as q(—x) = 
— q(x) and Aq(—x) = +Aq(x). 

• In the region — £ < x < £, one quark propagator represents a quark and the other one 
an antiquark. In this region, the OFPD's behave like a meson distribution amplitude. 

An alternative, but equivalent, parametrization of the matrix elements of the bilocal op- 
erators has been proposed by Radyushkin by expressing the longitudinal momentum 
fractions with respect to the initial nucleon momentum p instead of the average momentum 
P of Fig. 0. In this notation, the active quark has momentum k + = Xp + and A + = — (p + . 
The two variables X and ( are related to the previously introduced momentum fractions x 
and £ by: 

I = T3c7F- i = — 2 - (16) 

In terms of the variables X, ( and t, so-called 'non-forward parton distributions' (NFPD's) 
f q (X,(,t) were introduced by Radyushkin ||. They are equivalent to the OFPD's and are 
related to them by 



2 Note the typing errors in Eqs. (120, 121) of Ref. H, where the factors 2n should be replaced by 

4"7T. 
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H q (x^,t) = ^-j— {f q (X,(,t)} , for £<x<l, 
H q {x,i,t) = {f q (X,(,t) - f((-X,(,t)} , for 

H*(x,t,t) = ^ {-f*(C-X,C,t)} , for -1 <*<-£, (17) 

where / 9 are the quark and f q the anti-quark NFPD's respectively. 

The "double" nature of the skewed parton distributions according to whether < x < £ 
(0<X<C)or£<s<l(C<X< 1) was made more transparent by Radyushkin 
by expressing the quark momenta connected to the lower blob of Fig. |in a fraction x of 
the initial nucleon momentum p and a fraction y of the momentum transfer A as shown in 
Fig. |^. The relation with the previously introduced momentum fraction X is readily seen 
to be X = x + (y. In terms of the variables x , y, t, Radyushkin then introduced [pj|l4| the 
so-called double distributions F q , K q for quarks and F q , K q for the anti-quarks, according 
to: 



(p'|^(0) l + ^{z)\p) 

= N{p) 1 + N{p) ■ j^dx jf 1 dy {e~< £p+ - yA+ ) 2 "F 9 (£,|/,t) - e<* p+ + yA+ >~ F%x, y, £)} 



z+=z ± =0 

1 r l-x 



+ N(p')ia +K ^^N(p) 
2rriN 

-l-x 



x J dx dy \e- i ^ p+ - y ^ + ) z ~K q {i,y,t) - e< £p+ + ^ + ) z ~ K q {x, y, t)} , (18) 

with both x and y positive and x + y < 1. The double distributions are related to the 
NFPD's through integrals over y p|,|l4f : 



P(X, C, t) = dy F\X - (y, y, t) , for X > ( 

Jo 

rX/C 



f*(X,C,t)= dy F q (X-Cy,y,t), for X<C, (19) 

Jo 

with X = 1 — X and ( = 1 — (. Relations similar to Eq. ([I9| ) hold for the anti-quark distri- 
butions f q . Analogously to Eq. (|l"8D, one defines the helicity dependent double distributions 
F and K through the matrix element of a bifocal axial vector operator. 

Besides coinciding with the quark distributions at vanishing momentum transfer, the 
skewed parton distributions have interesting links with other nucleon structure quantities. 
The first moments of the OFPD's are related to the elastic form factors of the nucleon 
through model independent sum rules [0] . By integrating Eq. (|J) over x, one gets the 
following relations for one quark flavor : 

^ dx H«(x, e, t) = F q {t) , J + dx E\x, £, t) = F q {t) , (20) 
f % dx H*(x, £, t) = g\{t) , J + dx E«(x, £, t) = h\{t) . (21) 
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The elastic form factors for one quark flavor on the RHS of Eqs. (^(| [21]) have to be related 
to the physical ones. Restricting to the u, d and s quark flavors, the Dirac form factors are 
expressed as 

F u/p = 2F p + F n + pa ^ pd/p = 2 pn + pp + ps ^ ^ 

where Ff and F™ are the proton and neutron electromagnetic form factors respectively. The 
strange form factor is given by F° /p = F{, but since it is small and not so well known we set 
it to zero in the following numerical evaluations. Relations similar to Eq. (|22|) hold for the 
Pauli form factors F$. For the axial vector form factors one uses the isospin decomposition : 

gt = \gA + \&, gf = -\gA + \&. (23) 

The isovector axial form factor g& is known from experiment, with ^(0) ~ 1.267 ||15|| . For 
the unknown isoscalar axial form factor g^, we use the quark model relation : g\{t) = 
3/5 gA(t). For the pseudoscalar form factor h q A , we have relations similar to Eq. ([2~3|). 

The second moment of the OFPD's is relevant for the nucleon spin structure. It was 
shown in Ref. |3J that there exists a (color) gauge-invariant decomposition of the nucleon 
spin : 

l = J q + J 9 , (24) 

where J q and J g are respectively the total quark and gluon angular momentum. The second 
moment of the unpolarized OFPD's at t — gives 

1 f +1 

J q = - dxx [H g {x,£,t = 0) + E q {x,£,t = 0)] , (25) 

2 j_ x 

and this relation is independent of ^. The quark angular momentum J q decomposes as 

J q = ^AS + L q , (26) 

where AE/2 and L q are respectively the quark spin and orbital angular momentum. As AE 
is measured through polarized DIS experiments, a measurement of the sum rule of Eq. fl2"5]) 
in terms of the OFPD's, provides a model independent way to determine the quark orbital 
contribution to the nucleon spin and thus from Eq. (E3j) the gluon contribution. 

C. Modelization of the skewed parton distributions 

Ultimately one wants to extract the skewed parton distributions from the data but, in 
order to evaluate electroproduction observables, we need a first guess for them. We shall 
restrict our considerations to the near forward direction because this kinematical domain is 
the closest to inclusive DIS, which we want to use as a guide. Since E and E are always 
multiplied by the momentum transfer, this means that the amplitudes are dominated by H 
and H which we discuss first. 
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1. ^-independent ansatz for H and H 

The simplest ansatz, already used in our previous works is to write H and H as 

a product of a form factor and a quark distribution function, neglecting any £ dependence, 
which yields : 

H u ^(x, & t) = u(x) F 1 M/p (t)/2 , tf*/*(x, 1 1) = d{x) F? /P (t) , H s / p (x, £,t) = , (27) 

where m(x),c?(x) and s(x) are the unpolarized quark distributions. Obviously the ansatz of 
Eq. (|2"TD satisfies both Eq. (|T3| ) and the sum rule of Eq. (p0[), if one uses the valence quark 
distributions normalized as: 



-1 r+l 

dx uv(x) — I dx dy(x) = 1 . (28) 



o 



For all the calculations presented in this paper, we shall use the parametrization MRST98 
1 17] for the quark distributions. 

To model H q , we need the corresponding polarized quark distributions. For this we follow 
the recent work of Ref . [ IB| , where a next to leading order QCD analysis of inclusive polarized 
deep-inelastic lepton-nucleon scattering was performed and which yields an excellent fit of 
the world data. In this analysis, the input polarized densities (at a scale Ql = 1 GeV 2 ) are 
given by : 

Au v (x, Ql) = Vu A u x - 250 u v (x, Ql) , 
Ad v (x, Ql) = r] d A d x - 231 d v (x, Ql) , 

Aq(x, Ql) = n q A s x 5W S(x, Ql) , (29) 

where we assume a SU(3) symmetric sea, i.e. Aq = Au = Ad = As and S represents the 
total sea. On the rhs of Eqs. fl29D , the normalization factors A u ,Ad,As are determined so 
that the first moments of the polarized densities are given by i] u ,rjd,rjq respectively. For 
the valence quark densities, r\ u and rja were fixed in Ref. [0 by the octet hyperon (3 decay 
constants which yield : 

r+l r+1 

rj u = / dx Au v (x) w 0.918 , ij d = dx Ad v (x) w -0.339 . (30) 
Jo Jo 



The first moment of the polarized sea quark density was determined by the fit of Ref. flc 
which gives rjg w —0.054. 

Our ansatz for H q is to multiply the polarized quark distributions of Eq. (EH) by the 
axial form factor, neglecting any ^-dependence and sea contribution, which yields 

H«f*(x, t, t) = Au v (x) g u A /p (t)/g u A /p (0) , H d ^{x, £, t) = Ad v {x) g d J p (t)/g d J p (0) . (31) 

One can check that Eq. ( |T3] ) is verified by construction and that, using the (p3|) and the 
quark model relation to evaluate the axial form factors g U j[ p and g A p , the sum rule of Eq. (|2lD 
is satisfied within 10 %. 
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2. £- dependent ansatz for H and H 

To generate the dependence on £, we start from the double distributions introduced in 
Eq. fll8|) . As shown in Fig. |3|, the momentum flow in the double distributions is decomposed 
in a part along the initial nucleon momentum p and a part along the momentum transfer 
A . Therefore, a reasonable guess for the helicity independent double distribution F(x,y,t) 



was proposed in Ref. |TJ]], as a product of a quark distribution function q(x), describing 
the momentum flow along p, and an asymptotic "meson-like" distribution amplitude ~ 
— x — y), describing the momentum flow along A. Properly normalized, this yields : 



F*(x, y, t) = F?(t)/F?(0) q(x) 6 y \, y> , (32) 



y( l - g - y) 

[l-x) 

where similarly to Eq. (p7|), the t-dependent part is given by the form factor Fx. Similarly 
for the helicity dependent double distribution, this ansatz yields : 

F«(x, y, t) = g A (t)/g A (0) Aq(x) 6 ^'^'^ , (33) 



Starting from Eq.(^), the OFPD H q is then constructed by using Eqs. fllTD and (pT9|) . 
Analogously, the OFPD H q is constructed starting from Eq. (|33|). 

As an example, we show on Fig. |5| the dependence on £ predicted by this model for the 
valence (i-quark and total ci-quark skewed distributions. 



3. The skewed distributions E and E 

The modelization of E and E, which correspond to helicity flip amplitudes, is more 
difficult as we don't have the DIS constraint for the x-dependence in the forward limit. 
However, as already pointed out, E and E are multiplied by a momentum transfer and 
therefore their contribution is suppressed at small t. Nevertheless it was argued by the 



authors of Ref. |18] that the pion exchange, which contributes to the region — £ < x < £ of 
E, may be non negligible at small t due to the proximity of the pion pole at t = m\. We 



follow therefore the suggestion of Refs. |TH[T9| and evaluate E assuming it is entirely due to 



the pion pole. (By contrast we continue to neglect the contributions due to E since there 
is no dangerous pole in this case). According to this hypothesis, since the pion exchange is 
isovector, we must have 

fr/p = = I E ^ pole . (34) 

The t -dependence of E n ^ pole (x,!;,t) is fixed by the sum rule (pID pseudoscalar form factor 
h A (t) : 

+1 dxE^(x,£,t) = h A (t) — 9A(2M N f 
1 —t + m£ 

where E^ = E u / p — E d l p denotes the isovector OFPD and where we have approximated 
the pseudoscalar form factor by its pion pole dominance expression. 
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In the region — £ < x < £, the quark and antiquark couple to the pion field of the 
nucleon. Therefore, this coupling should be proportional to the pion distribution amplitude. 
For the later we take the asymptotic form, which is well supported now by experiment (see 



discussion in section |I V A|) . Expressing the quark's longitudinal momentum fraction z in 
the pion in the symmetric range — 1 < z < 1, the asymptotic distribution amplitude $ as 
is given by Q as {z) = 3/4 (1 — z 2 ), and is normalized as f\dz & a s(z) — !■ The light-cone 
momentum fractions of the quark and antiquark in the pion which couples to the lower blobs 
in Figs. ||, f|, are respectively given by (1 + x/£)/2 and (1 — x/£)/2. Therefore, E^_ po i e is 
finally modelled as 

K- pole = 9 (-£ < x < h A (t) i $ as f ~) , (36) 



which satisfies the sum rule of Eq. (|35|) . 

As there is no dangerous pole in the case of i£, it seems safe to neglect its contribution 
to the amplitude at small t, which of course does not mean that E itself is small. In this 
respect it is amusing to note that, if one adopts for E an ansatz similar to the one for H, 
Eq. (^), that is 

E u ^(x, £, t) = u(x) F 2 u/p (t)/2 , EV'(x, £, t) = d(x) F* /p (t) , E^ (x, £,t) = 0, (37) 

and if one uses this ansatz to evaluate the total spin carried by the quarks through the sum 
rule of Eq. (|25|) , one get the results shown in Table |. We can see that, even though there 
is little theoretical basis for the ansatz of Eq. fl3"?|) , we get the reasonable result that 43% of 
the nucleon spin in carried by the quarks. This is consistent with a recent QCD sum rule 
estimate |23[| (at a low scale) where the gluons were found to contribution to half of the 
nucleon spin. Note that only the information from the unpolarized parton distributions has 
been used to evaluate the quark contribution to the spin of the nucleon. 

By using the decomposition of the total quark spin J q as in Eq. ( |26| ) and by using the 



most recent values measured at SMC for the quark helicity distributions Am, Ad and 
As, one can extract the quark orbital angular momentum contribution to the nucleon spin. 
Using the factorized, ^-independent ansatz, the orbital contribution is estimated in Table | 
to yield about 15% of the nucleon spin. 



4- Other models of the skewed distributions 

Some model calculations of the skewed distributions exist. A bag model evaluation of 
the skewed distributions found (at a low scale) that they depend weakly on £. However, 
a calculation in the chiral quark soliton model |2l] found (also at a low scale) a strong 
dependence on £ and exhibits in particular fast "crossovers" at |x| = £. The chiral quark 
soliton model is based on the large- iV c picture of the nucleon as a heavy semiclassical system 
whose iV c valence quarks are bound by a self-consistent pion field. The crossovers in this 
model can be interpreted as being due to meson-exchange type contributions in the region 
— £ < x < £, which are not present in the bag model calculation and which are also not 
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contained in the phenomenological ansatz of Eqs. 



for the double distributions (see 



Ref. f2^] for a detailed discussion of these meson-exchange type contributions to the skewed 



quark distributions). It was noticed in Ref. [21 that this crossover behavior at 
lead to an enhancement of observables, which will be interesting to quantify. 



x 



£ could 



As a last remark, as we are mainly interested in this paper in giving estimates for 
electroproduction reactions in the valence region at only moderately large values of Q 2 
(1 ~ 20 GeV 2 ), which is where experiments of exclusive electroproduction reactions can be 
performed currently or in the near future, we neglect here the scale dependence of the skewed 
distributions. At very high values of Q 2 , the scale dependence should be considered. This 
scale dependence is described by generalized evolution equations that are under intensive 
investigation in the literature. 



III. LEADING ORDER AMPLITUDES FOR DVCS AND HARD 
ELECTROPRODUCTION OF MESONS 

A. DVCS 



Factorization proofs for DVCS have been given by several authors |2^,|2B|. They give 
the theoretical underpinning which allows to express the leading order DVCS amplitude in 
terms of OFPD's. Using the parametrization of Eq. (|j), the leading order DVCS tensor 
H'l.o.dvcs follows from the two handbag diagrams of Fig. |] as 



TT^U 

n L.O.DVCS 



1 r +1 

- \p»n v + p u n» - g^} / dx 

2 J -i 



+ 



x — £, + it x + £ — ie 



x 



12 DVCS 



(x,£,t) N(p)j.nN(p) + E p DVCS (x,Z,t) N(p)ia* xT ^N(p) 



+ \ l-ie^p K n x ] J + dx 



x — £, + ie x + £ — ie 



12 DVCS [ 



2,171 jv 



A • n 



x } £, t)N(p h.n l5 N(p) + E p DVCS (x, £, t)N(p ) l5 -—N(p) 



(38) 



with £ i23 = +1- We refer to Ref. for the formalism to calculate DVCS observables 
starting from the DVCS tensor of Eq. fl38|) . 

In the DVCS on the proton, the OFPD's enter in the combination 



9 9 9 



(39) 



and similarly for H, E and E. 

For the 7r° pole contribution to the DVCS amplitude, which contributes to the function 
E p DVCSl the convolution integral in Eq. ( j38|) can be worked out analytically. By using 



Eqs. (|3|, |36|) one obtains : 
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+1 



dx 



-i 



x 



IE 



- £ M*) 



(40) 



(|38|) , is exactly gauge invariant with respect 
However, electromagnetic gauge invariance 



The leading order DVCS amplitude of Eq. 
to the virtual photon, i.e. q u H£ u DVCS = 0. 
is violated by the real photon except in the forward direction. In fact q^H£ u DVCS ~ Aj_. 
This violation of gauge invariance is a higher twist effect of order 1/Q 2 compared to the 
leading order term H£ u DVCS . So in the limit Q 2 — > oo it is innocuous but for actual 
experiments it matters. Actually for any cross section estimate one needs to choose a gauge 
and this explicit gauge dependence for nonzero angles is unpleasant. In the absence of a 
dynamical gauge invariant higher twist calculation of the DVCS amplitude, we propose to 
restore gauge invariance in a heuristic way based on physical considerations. We propose to 
introduce the gauge invariant tensor Hp VCS : 



DVCS 



n L.O.DVCS 



[a ■ q 



r! U Xv 

Ix n L.O.DVCS 



(41) 



where is a four-vector specified below. Obviously Hp VCS respects gauge invariance for 
both the virtual and the real photon : 



DVCS 







% DVCS 







(42) 



Furthermore, as 



VA n L.O.DVCS 



(43) 



the gauge restoring term gives zero in the forward direction (Aj_ = 0) which is natural. We 
choose a? = p^ 1 because p ■ q is of order Q 2 , which gives automatically a gauge restoring 
term of order 0(1/Q 2 ). This choice for a M is furthermore motivated by the fact that in 
the derivation of the leading order amplitude of Eq. (|38|), only the p^ components at the 
electromagnetic vertices are retained 0. The above arguments lead to the following DVCS 
amplitude: 



TTfMU 

n DVCS 



H 



L.O.DVCS 



+ 



[p-q 



(A. 



H 



L.O.DVCS > 



(44) 



which is now gauge invariant with respect to both the virtual and real photons. We will 
illustrate the influence of this gauge invariance prescription, by comparing DVCS observables 
calculated with Eq. ( |44"|) and with the leading order formula of Eq. (^) . 



B. Hard electroproduction of mesons 

The OFPD's reflect the structure of the nucleon independently of the reaction which 
probes the nucleon. In this sense, they are universal quantities and can also be accessed, 
in different flavor combinations, through the hard exclusive electroproduction of mesons 
- 7r 0,± , rj, p 0,± , u>, 0, ... - (see Fig. ||) for which a QCD factorization proof was given in 



12 



Refs. ||f|]. According to Ref. ||], the factorization applies when the virtual photon is lon- 
gitudinally polarized because in this case, the end-point contributions in the meson wave 
function are power suppressed. Furthermore, it was shown that the cross section for a 
transversely polarized photon is suppressed by 1/Q 2 compared to a longitudinally polarized 
photon. 

Collins et al. || also showed that leading order PQCD predicts that the longitudinally po- 
larized vector meson channels (p ^, ujl, 4>l) are sensitive only to the unpolarized OFPD's 
(H and E) whereas the pseudo-scalar channels (7r 0,± ,^, ...) are sensitive only to the polar- 
ized OFPD's (H and E). It was shown in Ref. |T2| that the leading twist contribution to 
exclusive electroproduction of transversely polarized vectors mesons vanishes at all orders 
in perturbation theory. In comparison to meson electroproduction reactions, we recall that 
DVCS depends at the same time on both the unpolarized (H and E) and polarized (H and 
E) OFPD's. 

According to the above discussion, we give predictions only for the meson electropro- 
duction cross section by a longitudinal virtual photon. The longitudinal 72 + p — > M + p 
two-body cross section dai/dt is 



da L 
dt 



16vr 



m 2 N )A(s,-Q 2 ,m 2 N ) 2 



h N h' 



M 



0, h N ; h N ) 



(45) 



where h^, h M are the initial and final nucleon helicities and where the standard kinematic 



function A(x, y,z) is defined by 

A(x,y,z) = 

which gives A(s 



r 



y 2 + z 2 



2xy — 2xz — 2yz . 



(46) 



Q 



m 



N) 



2m N \q L \ 



where \qi\ is the virtual photon momentum in the 
lab system. As the way to extract da^jdt from the fivefold electroproduction cross section 
is a matter of convention, all results for dai/dt in this paper are given with the choice of 
the flux factor of Eq. (g§). 

In Eq. (^), M. L is the amplitude for the production of a meson with A m = by a 
longitudinal photon. In the valence region, the leading order amplitude is given by the 
hard scattering diagrams of Fig. ^ (Th part in Fig. |j). The amplitudes M. L for p° L and 
7T° electroproduction were calculated in Refs. [[J (see also Ref. | |271| ) for which the following 
gauge invariant expressions were found : 



4 1 






- — 




- / dx 


9 Q 


Jo z . 


2 J -i 



1 



x — £ + ie 



+ 



1 



x 



/€ 



H p pl (x, £, t)N(p'h.nN(p) + E^ (x, £, t)N(p )^ KA ^iV(p) 



(47) 



4 1 

— 26 - — 

9 Q 







- / dx 


Jo z . 


2 7-i 



+ 



X 



£ + ie x + £ 



ie 



A • n 



x(47ca s ) \ Hl (x,U)N(phM l5 N(p) + El (x,U)N(p) l5 —N(p)j, (48) 
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where & p {z) and $> w (z) are the p and 7r distribution amplitudes (DA) respectively. The 

), one 



factor 4/9 is a color factor corresponding with one-gluon exchange. From Eqs. (47, [ 
sees that the leading order longitudinal amplitude for meson electroproduction behaves as 
1/Q . As the phase space factor in the cross section Eq. ( |4~5|) behaves as l/Q 4 at fixed xb 
and fixed t, this leads to a 1/Q 6 behavior of da^jdt at large Q 2 and at fixed xb and fixed t. 

According to the considered reaction, the proton OFPD's enter in different combinations 
due to the charges and isospin factors. For the electroproduction of p°, which corresponds 
to the quark state (with spin S = 1) 



V2 

the corresponding skewed distribution is given by: 



p°) = 4={|uu> - \dd)} , (49) 



H :&^ - MI HU/P + l Hd/p ) ■ 



(50) 



For the longitudinal electroproduction of the u>l and 4>l vector mesons, the amplitudes have 
an expression analogous to Eq. (|47|) in terms of H% L and . By considering the uo and 
mesons as pure quark states (with S = 1) 

\u) = -^{M + 1^)} > \4>) = \ss) , (51) 
one then obtains for the corresponding OFPD's : 

HZfa&t) = ^{l HU/P ~ \ RdlP } > ( 52 ) 

Hl(x,U) = -\h^. (53) 

Besides the longitudinal electroproduction amplitude Eq. ( fi8|) for tt°, one again has anal- 
ogous expressions for other neutral pseudoscalar mesons such as the rj. Starting from the 
quark states (with spin S = 0) for pseudoscalar mesons and neglecting any mixing for the 
7] : 

k°) = -7={\uu) - \dd)} , (54) 
V2 

= -^={\uu) + \dd) - 2 \ss)\ , (55) 
V6 

one obtains for the polarized OFPD's : 

H P Ax,U) = ^ {j^ /p + \H d/ ^ , (56) 
H>(x,t,t) = ^= !^H U/P - \H d ^ + . (57) 

For the charged meson channels p and 7r ± , we obtain from \p + >= —\ud > and \p~ >= 
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\du > for the corresponding longitudinal electroproduction amplitudes (which have been 
given previously for the p in Ref. [28] and for the in Refs. p9|l9|) : 



M L n± 



4 1 

le 

9 Q 

i 



1 



dx 



x 



-2/3 



/V3 1 



-1/3 J x-^ + ie \ 2/3 J x + ^-ie_ 
(4vra s ) ^ 3 \x,^t)N(p) 1 MN(p) + E^(x,^t)N(p')ta Kxl ^N( P ) } , (58) 



— le 



4 1 
9 Q 



o 



fix 



-i 



z 

-2/3 \ 



i r 1/3 1 i 

-1/3 / x - £ + ie \ 2/3 Jx + ^-ie 

:(4vra s ) (^( 3 )(x,e,t)iV(p')7.n75^(p) + ^ (3 Hx,e,t)iV(p')75^iV(p) 
[ 2m N 



(59) 



where the w- and <i-quark charges appear in front of the direct and crossed terms and where 
the following isovector forms for the corresponding OFPD's enter : 



H (3 \x,Z,t) 
H (3) (x,Z,t) 



H d lv 
}J d /p 



(60) 
(61) 



In all of the above we have only given the expressions for the OFPD's H and H, but exactly 
the same isospin relations are valid for the OFPD's E and E. In particular, the isovector 
OFPD E^ provides a prominent contribution to the charged pion electroproduction am- 
plitude because it contains the t-channel pion pole contribution, given by Eq. fl3"6"|). As seen 
before for the ir° pole contribution in case of the DVCS, the convolution integral for the 
charged pion pole contribution in Eq. ( p9| ) can also be worked out analytically. In the case 
of the 7r + electroproduction amplitude, one has : 



+i 



dx 



+ 



-e d 



x 



£ + ie x + £ 
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E^p i e {x^,t) = — (e u -e d ) h A (t) 



(62) 



In the amplitudes for longitudinal meson electroproduction Eqs.(|58"|, |59"D , the meson dis- 
tribution amplitudes $ enter. For the pion, recent data |30| for the 7r°7*7 transition form 
factor up to Q 2 = 9 GeV 2 , which will be briefly discussed in section |IV A| , support the 
asymptotic form : 



= V2U 6z(l - 



(63) 



with f-x = 0.0924 GeV from the pion weak decay. With this asymptotic DA for the pion, the 
charged pion pole contribution to the amplitude M.^± can be worked out by using Eq. (|62|). 
where we insert the pion pole formula of Eq. (RH) for the induced pseudoscalar FF hA{t)- By 
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using the PCAC relation qa/ ' f n = 9-kNN / m Ni where g n NN is the nNN coupling constant, 
we finally obtain for the pion pole part of the amplitude Mf^± '■ 

M L n+ (tt+ - pole) = ieV^QF^ (Q 2 ) -g™^ N(p') l5 N(p) , (64) 

v — t + 

where F n represents the pion electromagnetic FF. The leading order pion pole amplitude 
is obtained by using in Eq. fl64|) the asymptotic pion FF F£ s , which is given by (see also 
section [IV B| ) : 



FT (Q 2 ) = • (65) 



For the r], the transition form factor has also been measured in Ref. ||30|| , which also 
supports an asymptotic shape for the rj distribution amplitude : 

%{z) = f v 6z(l - z) . (66) 



For the normalization, we adopt the value used in Ref. [|30|, which in the notation of Eq. ( |66[ 
is given by : f v m 0.138 GeV. 

For the vector mesons, no experimental determination of the DA exists besides its nor- 



malization. Both recent updated QCD sum rule analyses |3T| , p2| and a calculation in the 
instanton model of the QCD vacuum |33| favor a DA for the longitudinally polarized p 
meson that is rather close to its asymptotic form. The calculations in both models differ 
however in the deviations from the asymptotic form. In all calculations for vector meson 
electroproduction shown below, we will also use the asymptotic DA for the vector mesons 

0> v (z) = f v 6z(l - z) , (67) 

with f p ps 0.216 GeV, f u ps 0.195 GeV and ps 0.237 GeV determined from the electro- 
magnetic decay V — > e + e~ . 

When evaluating the leading order meson electroproduction amplitudes at relatively low 
scales, we use the value of the strong coupling which was also used in the QCD sum rule 
analysis of the vector meson distribution amplitudes of Ref. : a s (fi = 1 GeV) ~ 0.56. At 
high values of Q 2 , the running of the coupling has to be considered. However, the average 
virtuality of the exchanged gluon in the leading order meson electroproduction amplitudes 
can be considerably less than the external Q 2 , which is therefore not the "optimal" choice for 
the renormalization scale. In the next section, we will study the inclusion of the transverse 
momentum dependence in the considered hard scattering processes and will then be able to 
adapt the renormalization scale to the average gluon virtuality. The corresponding strong 
coupling constant will then be calculated within the convolution integral, through an IR 
finite expression, which reduces to the standard running coupling at larger scales. 



3 Because we want to use the same convention for the DA (Eq. (p7|)) for all vector mesons, we 
changed our definition of f p compared to our earlier work, i.e. f p of this work corresponds with 
y/2 f p used in Refs. f||. 
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IV. CORRECTIONS TO LEADING ORDER AMPLITUDE : INTRINSIC 
TRANSVERSE MOMENTUM DEPENDENCE AND SOFT OVERLAP 

CONTRIBUTION 



A systematic study of higher twist corrections to hard exclusive processes is beyond the 
scope of the present work. Our limited goal is to model here two important mechanisms 
which give rise to power corrections to the leading amplitude. 

First, we will consider the intrinsic (or primordial) transverse momentum dependence 
of hard exclusive processes. When one neglects the intrinsic transverse momentum of the 
active quark, the hard exclusive electroproduction amplitudes can be written as a one- 
dimensional convolution of a hard scattering operator and a non-perturbative soft quantity 
which depends only on the longitudinal momentum fractions of the quark. This neglect of 
the parton intrinsic transverse momentum is exact only up to corrections of order 1/Q 2 . 

Second, for the meson electroproduction reactions, which contain a one-gluon exchange, 
soft 'overlap' mechanisms can compete with the leading order amplitude at lower Q 2 values. 
In the case of meson form factors (FF) both corrections have been studied in details by several 
groups. It turns out that including these corrections allows a quantitative interpretation of 
the available data for the FF at Q 2 values down to a few GeV 2 . We will briefly review these 
corrections for the 7r°7*7 transition FF and the pion electromagnetic FF for Q 2 values in the 
range 1 ~ 20 GeV 2 . Our aim is to introduce the notations and to set the stage to calculate 
these corrections for the hard electroproduction. We will see that the DVCS is analogous to 
the 7r°7*7 FF because the leading order diagrams involve no gluon exchange. The leading 
order diagrams for meson electroproduction, which proceed through a one-gluon exchange 
mechanism, are analogous to those for the pion electromagnetic FF. 



A. Intrinsic transverse momentum dependence in tt 7*7 FF and in DVCS 



The calculation of the 7r°7*7 transition FF will serve as our starting point to calculate 
the corrections due to intrinsic transverse momentum dependence for the case of DVCS. 
The leading order contribution to the 7r°7*7 transition FF is given by the diagrams of Fig. [7] 
which involve a one- dimensional convolution integral over the quark's longitudinal momen- 
tum fraction z. The corrections to the leading order amplitude have been estimated by 
various groups (see Ref. [34] for a recent comparative discussion). We start here from the 
method used in in Refs. |35] which is based on the modified factorization approach of Refs. 
p6| and which we will extend to the case of hard electroproduction reactions. 

In this modified factorization approach, the amplitude corresponding to the two diagrams 
of Fig. [7] is expressed as a three-fold convolution integral over the quark's longitudinal mo- 
mentum fraction z and its relative transverse momentum fc^of the hard scattering operator 
T H and the pion light-cone wavefunction ^ n (z, k±) : 

F n o^(Q 2 ) = J dz J ^^(z,k ± )T H (z,k ± ;Q 2 ) ■ (68) 

The pion wavefunction ^ n (z, k±) represents the amplitude to find a pion in a valence qq state, 
where one quark has longitudinal momentum fraction z and relative transverse momentum 
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k±.The pion light-cone wavefunction is denned by the following bilocal quark matrix element 
at equal light-cone time (y + = 0) : 



V^zM) = \ Idy-e'W)*- I d 2 y ± e~ ik ^ (0\d(0)^u(y)\n + ( P7r )) 



, (69) 

y+=0 



where we have written the isospin structure for a positively charged pion moving with large 
momentum and where the normalization factor 1 / y/6 = 1/ \f2N c corresponds to the 
Brodsky-Lepage convention f37j. The hard scattering operator T# in Eq. (|68|) is calcu- 
lated from the diagrams of Fig. [7] by keeping the k± dependence in the intermediate quark 
propagators. Its evaluation yields : 

r „ ( ,, i;Q ^(^) 2 ^{_^ + _^_}, m 

where the factor l/(3\/2) is due to the squared quark charges in the 7r° state (Eq. 
When the kj_ dependence is neglected in Tjj the transverse momentum integral in Eq. 
acts only on ^^(z, k±) and one finds back the usual expression for the 7r°7*7 FF in terms of 
the pion distribution amplitude & n (z), which depends only on z : 



F^(Q 2 )^ I dz^T H {z-Q 2 ). (71) 



o 

The pion distribution amplitude $ n is defined as 

>w? d 2 k- 



|(jT 3 *Az,k ± ) = ^-j= $ v (z,fi F ) , (72) 

where we have indicated explicitely the factorization scale y,p. For \ip much larger than the 
average value of the quark transverse momentum, depends only weakly on fip and we 
will neglect this dependence in the following. Using the asymptotic DA Eq. ( |B13D for $ w , one 
obtains from Eq. ([71]) the well known leading order PQCD result: 



iW 7 (Q 2 ) Q — " ^ • (73) 



To investigate the corrections to the 7r°7*7 FF when Q 2 is not asymptotically large, one 
needs an ansatz for the non-perturbative wavefunction ^/ n (z,k±) in order to perform the 
integral of Eq. (RSI) over both z and k±. For the soft meson wavefunction, we follow Ref. 



35fl and adopt a gaussian form for the dependence on fcj 



».(,,£,) = ±& K exp -i >L >, (74) 
2V6 cr zz 2 a zz 



which satisfies Eq. (|72|) . The parameter a in Eq. (|74D is related to the average squared 
transverse momentum (k\) of the quarks in the meson. For the pion, this free parameter 
has been fixed in Ref. |38j through the axial anomaly which leads to the condition 
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1 AT 

dz*„{z,k ± = 0) = -y- . (75) 

Jit 



Using the asymptotic DA for $„. in the gaussian ansatz Eq. ([73) for the condition 
Eq. (^5)) yields 



<7 = 8tt 2 / 2 « 0.67 GeV 2 . (76) 

Furthermore using the gaussian ansatz for fy n with an asymptotic DA, one finds that the 
probability P q q of the valence quark Fock state of the meson, which is defined as 

is given by P^ = 2n 2 f 2 /a = 1/4. For the average squared transverse momentum defined by 

= 4- r dz i s^ii*(^)i 2 , (78) 



P qq Jo J 16vr 3 

one obtains an average quark transverse momentum of \J (k±) = a/ct/5 ~ 0.366 GeV. 

In Fig. ^ we show the Q 2 dependence of the 7r°7*7 FF evaluated with the asymptotic 
DA for the pion. The CLEO data 0] which extend to the highest Q 2 values (~ 9 GeV 2 ), 
clearly approach the leading order PQCD result of Eq. ([7JJ). The correction to the 1/Q 2 
scaling behavior in the range 1-10 GeV 2 is rather well described - as found in Refs. |35j - 
by the calculation including the intrinsic k± dependence and using the gaussian ansatz of 
Eq. ([74]) for the pion wavefunction. This kj_ dependence results in a reduction of the 7r°7*7 
FF compared to the leading order result when Q 2 decreases. At a value Q 2 ~ 3 GeV 2 , this 
reduction is about 20 %. 

We now study the intrinsic transverse momentum dependence for DVCS, which is a 
higher twist correction to the leading order DVCS amplitude of Eq. fl38|) . Our strategy is to 
evaluate the power corrections due to the intrinsic transverse momentum dependence only 
for the leading twist terms in Eq. fl38|) proportional to H, H, E and E. Guided by the leading 
order DVCS amplitude, we expect the contribution from these four OFPD's to be dominant 
and do not consider other higher twist OFPD's in this work. Starting from the handbag 
diagrams of Fig. § for DVCS, we evaluate the amplitude by keeping the transverse momentum 
dependence of the quarks propagating between the two electromagnetic vertices, as was 
discussed above for the 7r°7*7 FF. In analogy to Eq. (|68|) for the 7r°7*7 FF, the calculation 
of the DVCS amplitude leads to a three-fold convolution integral over the longitudinal and 
transverse components of k (see Fig. ^ for the definition). 

Let us first consider the generalization of Eq. (^ - representing the lower blob in Fig. |2| 
- when the momentum k in Fig. |2| has both a longitudinal component (k + = xP + ) and a 
transverse component (k±). We parametrize the bilocal quark operators matrix element at 
equal light-cone time (y + = 0) as 

P+ I dy-e< xP+ y f d 2 y ± e~^ (p'|^(-|) 7 +^(|)[p) 

= H%x, k±,£,t) N(p'h + N(p) + E q term , (79) 
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and 



d 2 y ± e- ik ^ ( J9 '|^(-|) 7 + 75 ^(|)b) 



P + J dy-e< xP+ y 
= H g (x,k ± ,£,t) N(p)-f + j 5 N(p) + EHerm 



■V- 



y+=0 



(80) 



where we have introduced fc^-dependent OFPD's H q (x,k±,C,,t),E q ,H q and E q . Remark 
that in writing down Eqs. (ff9j [80]) , we use a gauge as discussed in section |II A| , for which 
the gauge link between the quark fields is equal to one. 

Similarily to Eq. (|72|), where integrating the meson wavefunction over k± leads to the 
meson distribution amplitude, we find back the OFPD H(x,l;,t) introduced in Eq. (^) by 
integrating the k±- dependent distributions H(x, kj_,£,t), that is 



H'{x,£,t) 



d 2 kj 
~~8^ 



H q (x,k ± ,£,t) 



(81) 



and similarly for E, H and E. Indeed, one readily verifies that by integrating both sides of 
Eqs. (79, |8(j) over k± and by using Eq . flSll) , one finds back the leading twist parametrization 
of Eq. for the lower blob in Fig. |[ 

Next we must evaluate the hard scattering operator for the handbag diagrams of Fig. |2|. 
We keep the k± dependence only in the denominators of the (massless) quark propagators 
where it is most important. The quarks propagating between the electromagnetic vertices 
in Fig. have four- momenta (k + q — A/2) in the direct diagram and (k — q + A/2) in 
the crossed diagram respectively. Therefore, the transverse momenta in the propagators are 
(k± — Aj_/2) for the direct diagram and (k± + Aj_/2) for the crossed diagram. This leads 
to the following generalization of the leading order DVCS amplitude of Eq. (|38|): 



IpV + p v n» - g 



IJ.U] 





r d 2 k ± 


L dx j 




8tt 3 







H p ( 

n DVCS y 



X 



x, k±, £, t) N(p ) r y.nN(p) + E term 
1 



+ 



+ [-ie^Pnnx] 





r d 2 k ± 


C dx j 




8tt 3 



X 



x — + ie x + £ + — ie 
H p DVCS {xM,i,t) #(p') 7 .n 75 iV(p) + Eterm 
1 1 



x 



IE 



X + ^ 



16 



(82) 



with 



(83) 



One readily sees from Eqs. (^3, |83|) that by neglecting the transverse momenta in the quark 
denominators compared to Q 2 and by using Eq. fl8~l|) , one finds back the leading order DVCS 
amplitude of Eq. (|38|) . 
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In order to estimate the DVCS amplitude of Eq. (|82|), one needs to model the k±- 
dependent OFPD's. We will model here only the OFPD's H (x, k±, £, t) and H, as they give 
the dominant contribution to the DVCS amplitude in the near forward direction. Therefore 
we will neglect the contribution from the OFPD's E and E. 

To generate the fc^-dependent OFPD's, we start from the double distributions which 
we model in the near forward direction by a gaussian ansatz by analogy with what was 
done in Eq. (^) for the pion wavefunction. This yields for the helicity independent double 
distribution F q : 



F*{x,y,k ± ,t) = F?(t)/F?(0)q(x)6 



y(l - x - y) 
(l-xf 

4tt 2 



x 



ay(l-x -y) 



exp 



A-1 



2 ay(l -x-y) 



(84) 



For the helicity dependent double distribution F q , our ansatz yields 

y(l-x- y) 



F q (x,y,k ± ,t) = g q A (t)/g q A (0)Aq(x)6 



4tt 2 



1 



ay(l-x-y) 



exp 



2 ay(l-x-y) 



• (85) 



When integrating Eq. (|8~4]) over k±, one finds back the double distribution F q of Eq. 
i.e. Eq. ( [84"1) satisfies the normalization 



F q (x,y,t) 



d 2 k 



F q {x,y,k ± ,t) 



(86) 



In the same way, when integrating Eq. (|85D over k±, one finds back the helicity dependent 
double distribution F q of Eq. (|33|). 

The gaussian ansatz of Eqs. (|84] , ^5| ) depends upon one free parameter a which is related 
to the average squared transverse momentum of the quarks in the nucleon. We take here 
the same value as determined for the pion in Eq. ([76]). Having modelled the fc^- dependent 
double distributions, the fc^-dependent OFPD's H and H are then obtained from Eq. (53) 
using the k±_ dependent analogues of Eqs. ([H]) and ([19]). 

When calculating the k±- dependence for the hard meson electroproduction reactions in 
section [IV B| , we will find it convenient to introduce the OFPD's in impact parameter space. 
They are obtained by taking the Fourier transform of the dependent OFPD's with respect 
to k : 



H q (x,b,£,t) 



d 2 k ± 
(2^ 



(87) 



where b plays the role of a transverse distance between the active quarks. Using Eq. (|ST|). 
one readily verifies that the ordinary OFPD's correspond to a zero transverse distance : 
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7 }-H'(x,b = 0,Z,t) = H*(x,£,t)- (88) 

We show in Fig. |9| the impact parameter dependence for the valence down quark OFPD H v 
using the ansatz of Eq. (|84]) . Note that the values at the origin {b = 0) in Fig. || can be 
directly read off from the corresponding ordinary OFPD displayed in Fig. H. 



B. Intrinsic transverse momentum dependence in pion electromagnetic FF and in 

hard meson electroproduction reactions 

Having discussed the intrinsic transverse momentum dependence of the DVCS amplitude, 
we now turn to the hard meson electroproduction reactions. As before, where we have made 
the parallel between the 7r°7*7 FF and DVCS, we will find it very instructive to study the 
intrinsic transverse momentum dependence of the hard meson electroproduction amplitudes 
by comparison with the charged pion electromagnetic FF. In contrast to the 7r°7*7 FF, 
the leading order contribution to the pion electromagnetic FF proceeds through one-gluon 
exchange as shown in Figs. [XOj a, b. 

To calculate the form factor, we consider the process in a frame (e.g. Breit frame) where 
the initial (final) pion moves with large momentum along the positive (negative) z-direction 
and where the momentum transfered by the virtual photon points along the negative in- 
direction. We denote the longitudinal momentum fractions of the quarks in the initial 
(final) pion by x (y) and their relative transverse momenta by k± (l±). The expression of 
the pion electromagnetic form factor including the transverse momentum dependencies is 
given by the six-fold convolution integral : 

FAQ 2 ) = f dy j j\ x j ^*AyJ±)T H (yJx,xM;Q 2 )*Ax,k x ) . (89) 



The hard scattering operator in Eq. (p9[) can be easily calculated by the diagrams of Figs. [L0 



a, b and two other diagrams where the virtual photon couples to the lower quark line. It 



has been shown in Ref. |4(J that the most important transverse momentum dependence 
in Tjj comes from the gluon propagator. The reason is that if one neglects the transverse 
momentum dependence, the gluon virtuality then depends quadratically on the quark lon- 
gitudinal momentum fractions (x and y). In contrast, the virtuality of intermediate quark 
lines depends only linearly upon the longitudinal momentum fractions. Therefore the gluon 
propagator weights most heavily the end-point regions and it is important in these regions 
to include the transverse momentum dependence to get a consistent PQCD calculation. We 
keep therefore the transverse momentum dependence of T# only in the gluon propagator, 
which simplifies considerably the integral of Eq. ( |55| ) . Using the symmetry of the integrand 
in Eq. (|89|) under x <-> x and y <-> y, this yields for T# : 

THiy , r ± , x, s ±i q 2 ) - I (4™,) xyQl + \ h -t 1? • < 90 > 

where the factor (4/9) is a color factor due to one-gluon exchange. Neglecting the transverse 
momenta dependence in T# and using asymptotic DA's, one obtains the well known leading 
order PQCD result for F n : 
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Q 2 



(91) 



The calculation of the full six-fold convolution integral of Eq. (p9|) can be simplified by 
noting that the hard scattering operator Tjj of Eq. (^) depends only on the transverse 
momentum difference m± = k± — Therefore, one has advantage to express Eq. fl39l) in 
impact parameter space. The pion wavefunction in impact parameter space is obtained from 
ty n (x, k±) by the Fourier transform 



d 2 k ± 



D i k±-b 



(92) 



where b represents the transverse separation between the quarks in the pion. For the gaussian 
ansatz of Eq. (0), ^f n (x, b) is given by : 



^(X, b) = 4:71 



(x 



2v/6 



exp 



- - ab 2 xx 
2 



(93) 



The hard scattering operator in impact parameter space is given by 

d 2 mj_ 



T H (x,y,b;Q 



(2t) s 



e^ b T H (x,y,rh ± ;Q 2 ) 



4 i 

- (l6na a (fi 2 R )) — Ko(y/xybQ) , 
3 Zix 



(94) 



where the strong coupling is evaluated at the renormalization scale fin to be specified shortly, 
and where the modified Bessel function K$ in Eq. ( |94"D is obtained by Fourier transforming 
the gluon propagator (in Eq. (^)) with spacelike four-momentum. As and T# depend 
only on the magnitude (b) of b, Eq. (|89| ) for the pion FF reduces to a three-fold convolution 
integral : 



i f i 

dx / dy 
ii Jo Jo 



dbb 



S 7T 



^{x,b)T H {x,y,b; Q 2 ) V n (y,b) . 



(95) 



To calculate the hard scattering operator of Eq. (|94"D, one has to specify the renormalization 
scale /i/j in the strong coupling constant. Here we follow the recent work of Ref. [^l| and 
take it as n 2 R = ^Q 2 , which represents in a sense the average gluon virtuality in the leading 
order diagrams for the pion FF. Near the endpoints, where the quark longitudinal momenta 
vanish (x ~ or y 0), the gluon virtuality becomes very small and one leaves the region 
where the PQCD result for the running coupling can be used. An often used practise is to 
freeze the coupling at low gluon virtuality at some finite value (~ 0.5). Instead of this ad- 
hoc procedure, the pion FF was evaluated in Ref. |5T| using the infrared (IR) finite analytic 
model for a s of Refs. ||2,[13]. In Ref. |[42|| , Shirkov and Solovtsov have deduced from the 



asymptotic freedom expression for a s and by imposing Q 2 analyticity, the following IR finite 
result for the (one-loop) strong coupling : 



a s (Q ) = ir 
Po 



^Q 2 /A 2 QCD 



+ 



A 2 



A 2 



Q 1 



(96) 
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where (3q = 11 — 2/3 and Aqcd is the QCD scale parameter. The first term in Eq. fl96|) 
is the standard (one-loop) asymptotic freedom expression and the second term ensures that 
the coupling has no ghost pole at Q 2 = Aq CD . At Q 2 = 0, the coupling constant takes on 
the finite value ot s {Q 2 = 0) = 47r//3 (~ 1-40 for Nj = 3), which depends only upon group 
symmetry factors. Therefore, Eq. fl96|) provides a coupling that can be used from low to high 
values of Q 2 without any adjustable parameters other than Aqcd- The value of Aqcd in the 
one-loop case can be fixed from a s (M 2 ) « 0.34 and results in Aqq d « 280 MeV - where 
the superscript refers to the analytic model of Eq. fl96|) (it was noted in Ref. [42] that this 



value of Aq CD corresponds with Aqq d ~ 230 MeV when using the MS scheme in the usual 
renormalization group solution). As an example one obtains from Eq. fl96|) a°, n {l GeV 2 ) ~ 
0.43. Having fixed Aqcd, the convolution integral of Eq. (P5|) for the pion form factor can 
now be evaluated with the IR finite strong coupling of Eq. ((5(J), as proposed recently in Ref. 
0- 

The result is shown in Fig. [11], where the leading order PQCD expression of Eq. ( pl|) 
for the pion form factor (dotted line) is compared with the result including the intrinsic 



transverse momentum dependence (dashed-dotted line). As is seen from Fig. 11, the leading 
order PQCD result is approached only at very large Q 2 . The correction including the 
transverse momentum dependence gives a substantial suppression at lower Q 2 (about a 
factor of two around Q 2 ~ 5 GeV 2 ). At these lower Q 2 values, the inclusion of the transverse 
momentum dependence renders the PQCD calculation internally consistent in the sense that 
the dominant contributions come from regions in which the coupling is relatively small. In 
addition to the reduction due to the transverse momentum dependence there is an additional 
suppression due to the Sudakov effect which ensures that at large Q 2 , soft gluon exchange 
between quarks is suppressed due to gluonic radiative corrections. At the intermediate 
values of Q 2 shown in Fig. [11], the exponential reduction due to the Sudakov form factor was 
calculated and found to be small provided one has already taken into account the intrinsic 
transverse momentum dependence, as was also noted in Ref. Only at very large values 
of Q 2 , the Sudakov form factor takes over and yields an additional reduction compared with 
the one due to the intrinsic transverse momentum dependence. As we give in this paper 
only predictions for low and intermediate Q 2 values, we will not consider the Sudakov effect. 

We now transpose the result for the pion electromagnetic FF and study the transverse 
momentum dependence in hard meson electroproduction amplitudes. Keeping the transverse 
momentum dependence of the hard scattering operator only in the gluon propagators, we 
obtain for the longitudinal electroproduction amplitude of p° L : 
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When neglecting the transverse momenta dependencies k± and l± in the gluon propagators 
in Eq. (|97|) and by using the normalization conditions for the p wavefunction (analogous as 
Eq. (|72| ) for the pion) and the normalization condition for the OFPD's (Eq. fl3T|)), we find 
back the leading order p\ electroproduction amplitude of Eq. (|47|). 

The evaluation of Eq. (|97|) requires to perform a six-fold convolution integral in analogy 
with Eq. (^) for the pion electromagnetic FF. To make the calculation tractable, we will 
limit ourselves here to the forward direction where Aj_ = 0. In this case, the gluon prop- 
agators depend only on the transverse momentum difference k± — l± and one can use the 
impact parameter space, as for the pion form factor. Using Eq. (^) for the wavefunction 
and Eq. (^) for the OFPD, Eq. (|97|) for the hard p\ electroproduction amplitude becomes 
in the forward direction (Aj_ = 0) : 



M L pl (k A 



4 1 (Q 
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H p p(L {x,b,U)N{p)j.nN{p) + EtermJ (4vra s (/4)) - T H (x, £, z, 6; Q 2 ) , (98) 
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where T H (x,£,, z,b;Q 2 ) represents the hard gluon propagators in impact parameter space. 
It is given by 

T H (x,£,z,b; Q 2 ) 
1 f°° 

= — / dm±m±Jo(bm±) 

= ^ if W (yz (x - 0/(20 Q bj + Ko(yz(x + Z)/{2()Qb)} , for x > £ , (99) 

= ^ {-K (y Z ^-x)/{20Qbj + K (yz{x + Z)/(2QQb)} , for < x < £ , (100) 

where the analytical results are obtained through a contour integration in the complex plane. 
The modified Bessel function K originates from Fourier transforming a gluon propagator 
with spacelike momentum which appears in the direct diagram when < x < £ and in the 
crossed diagram. The Hankel function in Eq. (p9|) originates from Fourier transforming 
a gluon propagator with timelike momentum, which appears in the direct diagram when 
x > £. In Eqs. fl9"9| , |100| ), the analytical expressions of T# are given for positive values of x. 
The expression of Tjj for negative values of x is directly seen to be minus the expression for 
positive values of x. 

For the estimates of meson electroproduction, which include the transverse momentum 
dependence, we use the IR finite strong coupling of Eq. (|96[) , as in the case of the pion form 
factor. The scale pn at which the strong coupling a s in Eq. (p8|) is evaluated, is given by 
the largest mass scale in the hard scattering, which we can take as 



max 



„ Q 2 i 



'101) 



for the direct diagram and 
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H% = max[\x + t\z^, ^) , (102) 



for the crossed diagram. 



C. Soft overlap mechanism for pion electromagnetic FF and for meson 

electroproduction reactions 

For hard exclusive reactions that proceed through one-gluon exchange at leading order, 
soft 'overlap' mechanisms can compete with the leading order amplitude at intermediate 
values of Q 2 . We will estimate the soft overlap contribution to meson electroproduction by 
analogy with the soft overlap to the pion electromagnetic FF. 

The overlap contribution to the pion electromagnetic FF is shown in Fig. pTD|c. Its cal- 
culation is most easily performed in the frame where the spacelike virtual photon has only 
a transversal momentum q± (i.e. q + = and Q 2 = qf_). Although the final result for the 
FF cannot depend on the reference frame in which one performs the calculations, the cal- 
culation however in a frame where q + 7^ is more complicated due to additional Z-graph 
contributions (see e.g. Ref. |46| for a comprehensive discussion). 

Indeed, in a frame where q + 7^ 0, one obtains two types of contributions to the FF. The 
first corresponds to the light-cone time ordering in the diagram of Fig. |l0| c, where the initial 
pion consists of partons with light-cone momenta k + and (pjjr — k + ) before the interaction, 
and where the final pion consists of an active parton with momentum (k + + q + ) and a spec- 
tator parton of momentum (p+ — k + ). In this case, the non-perturbative objects at both 
meson sides are the pion light-cone wavefunctions *f? n as defined before. 

The second type of contribution corresponds to the situation where the virtual photon 
creates a quark- ant iquark pair with light-cone momenta (q + — k + ) and k + prior (in light-cone 
time) to the interaction with the bound state. In contrast to the first type of contribution, 
these so-called Z-graph contribution cannot be expressed as the product of two pion light- 
cone wavefunctions. Indeed, it involves the amplitude for a parton to split into a meson 
and another parton, which is a different non-perturbative object than a meson light-cone 
wavefunction. However, as pair creation from a photon with zero momentum q + = is not 
possible for light-cone time ordered diagrams, these pair creation or Z-graph contributions 
vanish when q + = 0. Therefore, it is very convenient to estimate the soft overlap contribu- 
tion to the pion FF in the frame q + = where only the first contribution survives, and which 
is given by the Drell-Yan formula as a product of two meson light-cone wavefunctions : 



F^\Q 2 ) = J dz j ^y n ( z ,k ± + zq) *„(z,k ± ) . (103) 

In the arguments of the soft meson wavefunctions in Eq. ( |103|) , the relative transverse 
momenta enter. Using the gaussian ansatz of Eq. ( [74] ) for the soft pion wave function ^ in 
Eq. ( |103j ), one obtains : 
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The soft overlap contribution of Eq. (|104|) to the pion electromagnetic FF is shown in Fig. [11], 



calculated with the asymptotic distribution amplitude for It is seen from Fig. that 
although the soft overlap contribution drops faster than 1/Q 2 , one has to go to rather large 
values of Q 2 before the perturbative one-gluon exchange contribution dominates. In the 
intermediate Q 2 range shown in Fig. O, the soft overlap contribution makes up more than 



half of the total form factor. One also sees from Fig. [Il| that the sum of the soft overlap 
and the perturbative one-gluon exchange contributions, corrected for transverse momentum 
dependence, is compatible with the existing data. When comparing to the reported data for 
the pion electromagnetic FF, one has to be aware however that these sparse data are rather 
imprecise at the larger Q 2 values, as they are obtained indirectly from pion electroproduc- 
tion experiments through a (model-dependent) extrapolation to the pion pole. 

We now calculate the soft overlap contribution to longitudinal electroproduction of 
mesons, which is shown in Fig. [12]. As we are interested in this paper in the study of 
meson electroproduction reactions at a non-zero transfer A, the choice of frame to suppress 
Z-graph contributions completely is not as obvious as in the form factor case. Among the 
natural choices, one has now the possibility to choose a frame where q + = or a frame 
where A + = 0. As we are concerned in this work with the kinematical regime of large Q 2 
and small momentum transfer t = A 2 , it seems natural to estimate the overlap contribution 
in a frame where q + = 0, which will suppress Z-graph contributions at the virtual photon 
vertex. Remark that the choice of frame q + = is different from the one used before in the 
calculation of the leading order amplitude for hard meson electroproduction. We will use 
this different choice of frame only to provide an estimate for the soft overlap contribution 
and will therefore show the results for the leading order and overlap amplitudes separately. 
To add the leading order and overlap amplitudes coherently, one would have to perform a 
Lorentz boost. 



Making the choice q + = 0, we show in Fig. \LQ the leading (+, _L) components of the 
external and quark momenta in the overlap diagram, where the + components refer here 
to the initial nucleon momentum, to keep a close analogy to the overlap calculation for 
the pion FF. In the present work, we will only evaluate the contribution from the region 
< x < C in Fig. 12, in which case the meson vertex can be parametrized by the usual 
meson wavefunction. In order to estimate also the contribution outside this region, when 
x > (, we would again need to know how to parametrize the nonperturbative object where 
a quark splits into a quark and a meson. We postpone to future work the investigation of 
this problem. 

In the region < x < ( we can parametrize the valence quark state of a longitudinally 
polarised vector meson as: 

m = tJ? z I tStS (*■ Sj 71 1 * * + * * > ' (105) 



where is a color normalization factor and and q^ are the quark and antiquark 

states in Fig. [12| with helicities A = +1/2 (|) and A = —1/2 (J,). The quark states are 



normalized as |37j 



(q x ,(k' + ,k' A _)\q x (k + ,k ± )) = 2k + (2n)H(k' + -k + ) 5 2 {k' L -k^j S yx . (106) 
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The light-cone valence wavefunction tyy L in Eq. (|105|) , depends upon the quark relative 
light-cone momentum fraction (z) and the quark relative transverse momentum (k±) in the 
meson. This valence wavefunction tyy L is given, in analogy with Eq. (|69|) for the pion, and 
for a meson with py± = 0, by the bilocal quark matrix element at equal light-cone time 
(y + = 0): 



*v L (zj±) = 4 J dy-ei z tiy J d 2 y ± e-^ (0|g(0) 7 + q(v)\V L (p v )) 



• (107) 

y+=0 



The flavor structure for different vector mesons (p° , , to , 0) is understood implicitely in 
Eqs. (TIU51, |TO7D. 

We next have to evaluate the lower blob in the overlap diagram of Fig. [1^. As we 
restrict ourselves to the contribution from the region < x < (, this lower blob corresponds 
to the amplitude to find a qq state at the nucleon side. In the calculation of the leading 
order meson electroproduction amplitude, when using a frame where the initial nucleon's 
transversal momentum is zero, this amplitude is parametrized by the 'meson-like' part of the 
NFPD's. In the q + = frame with non-zero transversal momentum q±, the relative quark 
transverse momentum will appear as an argument in the fc^-dependent NFPD. Therefore, 
in the < x < ( region the lower blob in Fig. is : 

^ , f d 2 k\ „, - x _ . 1 

2C q±,U) l\f x (C-x) p n 



N(p )l + N(p) ■ ^ dx I —±f( x ,k ± + —q ± ,C,t) 7 — ^ = , \qm + gifr) • (1° 8 ) 



In Eq. ( |108|) , / is the k±- dependent unpolarized NFPD (as we are interested to calculate 
longitudinally polarized vector meson electroproduction), which is related to the ordinary 
NFPD as in Eq. flS!]). The term corresponding with the OFPD E is neglected here. Remark 
that in the infinite momentum limit, where one neglects the quark transverse momentum 
compared with its longitudinal momentum, Eq. ( |108| ) yields the Lorentz structure ~ (7") 
that appears in the leading order amplitude of Eq. 

We can now calculate the amplitude of the soft overlap mechanism by taking the matrix 
element of the electromagnetic current J M between the vector meson state \ Vl) and the \qq) 
state at the nucleon side of Eq. (|108[). This yields : 



My L , soft = (-te)N(ph + N(p) 

ft f d 2 k± x 1 1 

x / dx / TT f( x ' k± + ^7^' *) ~A ~r=n \ T < Vl I J ' e I ^1 + > ' ( 109 ) 
Jo J 2Q 4 yj x (C - x) p + 

where e is the photon polarization vector. By using the vector meson state of Eq. (|l05l) , 
Eq. ( |109| ) can be worked out as 

M Vl ,soft= He)iV(p') 7 + iV(p) 



x r £ 



* I dx J f(x,k ± + —q ± X,t)^y VL \z=^k ± + (l-^)q ± 







x <^ ( q^(xp + , k ± + q ± ) I J • e | g T (xp + , fcjj) 



+ (^(^ + ,^ + ?l)U-e|?;(^p + ,^±))] , (HO) 
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where the relative quark momentum arguments enter in the vector meson valence wave- 
function and in the NFPD. As we only give predictions in this work for the longitudinal 
electroproduction amplitude, we only keep the good current component J + , i.e. J • ejj — 
J + ■ (where is the polarization vector for a longitudinal photon). For the current 
operator J + , the quark matrix elements in Eq. ( |110| ) yield : 

( qy{xp + , k± + q ± ) | J + | q x (xp + , k±)) 
= u(xp + , kj_ + q±_] A') 7 + u(xp + , kj_; A) = (2xp + ) 5 X y ■ (111) 
With Eq. Ql 11 ), the soft overlap amplitude for production of a longitudinally polarized vector 



meson by a longitudinal photon can finally be written as : 

Mb L , sof t= He) e~ L N{p'h + N{p) 

x-L^dx J ^f {x7 k ± + ^ ±X: t)^ VL ^^J ± + (l-^)q^ . (112) 

In the frame q + = considered here, g M = (0,q±,0) and = (1,0,0,0), which yields 
e L = l/V2. 

To evaluate the soft overlap formula Eq. ( |1 12[ ) for meson electroproduction, we need to 



model the fc^-dependent NFPD for x < (. We do that by analogy with the gaussian ansatz 
of Eq. ( |7^ ) for the soft part of the meson wavefunction, that is: 



f(x,k' ± ,(,t) = f(x,(,t) % v ex P"I T } ' X ^C, (H3) 



where uat is now related to the average squared transverse momentum of the quarks in 
the meson state at the nucleon side. The ordinary unpolarized NFPD f(x,(,t) entering 
in Eq. ( |113| ) is modelled as discussed before through its definition Eq. (|T9"D in terms of the 



double distribution and through the ansatz of Eq. ( p2|) for the double distributions. With 
the ansatz of Eq. ( |113| ) for the NFPD and Eq. ([F3]) for the soft meson wavefunction (with 
parameter <jy L for the longitudinally polarized vector meson), the integral over transverse 
momentum in the soft overlap formula of Eq. ( |112| ) can be worked out analytically as : 

M l Vl >soft = He) -L iV>') 7 +iV(p) 

*fv L ^— f'dxFixXit) exp(- Q2 , {1 ~l /0 ), (H4) 



a N + a VL Jo { 2(a N + a VL ) x/Q 

where fy L appears in the asymptotic meson distribution amplitude as given before in 
Eq. d57|) . Remark that Eq. ( [L14|) is a generalization of the overlap formula of Eq. (|104| ) for 
the pion FF using a gaussian ansatz for the transverse momentum dependence. Note also 
that the overlap amplitude is purely real. In contrast, the leading order meson electropro- 
duction amplitude, including corrections for the intrinsic transverse momentum dependence, 
is complex and is furthermore dominated by its imaginary part. Therefore, one can expect 
to find qualitative differences with respect to the case of the form factor. 

In the actual calculations of this work, we use as first guess for oy h and a^, the same 
value as found before (Eq. ([76])) for the pion : o Vl w a N w 0.67 GeV 2 . In principle, the 
parameter could be fixed independently however, if data are available. 
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V. RESULTS AND DISCUSSION 



In this section we present results for several observables for meson electroproduction and 
DVCS. We will show the leading order predictions and the effects of the power corrections 
discussed in the previous section. In our previous work we gave predictions for the 

leading order amplitudes using a simple ^-independent factorized ansatz for the OFPD's. 
Therefore, the ^-dependent ansatz used here allows to study quantitatively the dependence 
on the skewedness in the OFPD's. 



We first show in Fig. [13] the angular and energy dependence of the leading order me- 
son electroproduction and DVCS differential cross sections in kinematics accessible at JLab, 
HERMES and COMPASS. For the electroproduction of photons, one cannot disentangle the 
DVCS from the Bethe-Heitler (BH) process where the photon is radiated from an electron 
line. In our calculations we add coherently the BH and DVCS amplitudes. From a phe- 
nomenological point of view, it is clear that the best situation to study DVCS occurs when 
the BH process is negligible. For fixed Q 2 and xb, the only way to favor the DVCS over 
the BH is to increase the virtual photon flux and this amounts to increase the beam energy. 
This is seen on Fig. [13] where we compare separately the BH, the DVCS and the coherent 
cross section for different beam energies. According to our estimate, the unpolarized (I, I 7) 
cross section in the forward region is dominated by the BH process in the few GeV region. 
To get a clear dominance of the DVCS process one needs a beam energy in the 100 GeV 
range. In the valence region (xb ~ 0.3) and for a Q 2 value of Q 2 ~ 2.5 GeV 2 , one sees that 
at = 200 GeV, the DVCS cross section is about two orders of magnitude larger than the 
BH in the forward direction. By going to smaller xb, the BH cross section increases. This 
is due to the fact that in the BH process the exchanged photon has 4-momentum (q — q') 
which gives al/t behaviour to the amplitude. The value of t in the forward direction (t min ) 
becomes very small for small x b values. The resulting sharp rise of the BH process in the 
forward direction at small xb puts therefore a limit on the region where the DVCS can be 
studied experimentally. At COMPASS kinematics, xb — 0.1 seems to be the lower limit. 
Although the BH is not a limiting factor at high xb, one cannot go too close to xb = 1 in 
order to stay well above the resonance region. 

In Fig. 13|, the leading order predictions for the p° an d tt + fivefold differential electropro- 
duction cross sections are also compared in the same kinematics and one again observes that 
the virtual photon flux boosts the cross sections as one goes to higher beam energies. In the 
valence region (i^ ~ 0.3), the p° L cross section, is about an order of magnitude larger than 
the DVCS cross section. For the 7r + electroproduction cross section, one remarks the promi- 
nent contribution of the charged pion pole (OFPD E), which gives a flatter t-dependence 
than the contribution of the OFPD H, because the E contribution comes with a momentum 
transfer in the amplitude. 

In Fig. [13] we compare in COMPASS kinematics, the leading order predictions for p° L 
electroproduction using the ^-dependent ansatz and a ^-independent ansatz which we used 
in our previous work MM. Both calculations in Fig. [T3| use the MRST98 parton distribu- 
tions |HJ as input and thus the OFPD's in both cases reduce to the same quark distributions 
in the forward limit and satisfy the first sum rule. One sees that the cross sections differ 
by a factor of two, which indicates that the skewedness of the OFPD's a priori cannot be 
neglected in the analysis. 
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Although it is clear from Fig. [13] that a high energy beam such as planned at COMPASS 
is preferable, one can try to undertake a preliminary study of the hard electroproduction 
reactions using the existing facilities such as HERMES or JLab, despite their lower energy. 
Concerning JLab, the high luminosity available compensates for the low cross section, and 
its good energy and angular resolution permits to identify in a clean way exclusive reactions. 
To make a preliminary exploration of the reaction mechanisms in the few GeV regions and 
to test the onset of the scaling, the measurement of the p° L leptoproduction through its decay 
into charged pions seems the easiest from the experimental point of view as the count rates 
are the highest. An experiment to explore the pi electroproduction at JLab at 6 GeV with 
the CLAS detector has been proposed 0. 

For the 7 leptoproduction at low energies, we suggest that an exploration of DVCS might 
be possible if the beam is polarized. The electron single spin asymmetry (SSA) does not 
vanish out of plane due to the interference between the purely real BH process and the 
imaginary part of the DVCS amplitude. Therefore, even if the cross section is dominated 
by the BH process, the SSA is linear in the OFPD's. To illustrate the point, we show on 
Fig. [15] the unpolarized cross section for a 6 GeV beam and the SSA at an azimuthal angle 
<f) = 120°. When the angle between the real and virtual photons is in the 0° - 5° region, a 
rather large asymmetry is predicted, even though the cross section is dominated by the BH 
process. We also show on Fig. [15] the effect of the ir° pole in the OFPD E to the DVCS. One 
sees that at small angles (small t), the effect of the 7r° pole is quite modest and increases at 
larger angles. 

We furthermore illustrate the effect of the gauge restoring term for the DVCS amplitude 
in Eq. (fHj). We have plotted in Fig. [15] the SSA both for the gauge invariant and non gauge 
invariant amplitudes. For the non gauge invariant amplitude of Eq. ([38]), the SSA is shown 
both in the radiative gauge and in the Feynman gauge. As expected, all predictions are 
identical at small angle. At larger angles the gauge dependence clearly shows up, especially 
in the Feynman gauge. 

As the SSA accesses the imaginary part of the DVCS amplitude, the real part of the 
BH-VCS interference can be accessed by reversing the charge of the lepton beam since this 
changes the relative sign of the BH and DVCS amplitudes. We have given in Ref. || an 
estimate for the e + e~ asymmetry at 27 GeV, which yields a comfortable asymmetry in the 
small angle region. This may offer an interesting opportunity for HERMES, although the 
experimental (e + e~) subtraction might be delicate to perform. 

Before considering the extraction of the OFPD's from electroproduction data, it is com- 
pulsory to demonstrate that the scaling regime has been reached. In Fig. [16], we show the 
forward longitudinal electroproduction cross sections as a function of Q 2 and compare the 
L.O. predictions for different mesons. The leading order amplitude for longitudinal electro- 
production of mesons was seen to behave as 1/Q. Therefore, the leading order longitudinal 
cross section da^/dt for meson electroproduction behaves as 1/Q 6 - In Fig. [16], we give all 
predictions with a coupling constant frozen at a scale 1 GeV 2 as explained in section PVB . 
When using a running coupling evaluated at the scale Q 2 , one probably underestimates the 
cross section in the range Q 2 ~ 1 - 10 GeV 2 , as the average gluon virtuality in the L.O. 
meson electroproduction amplitudes is considerably less than Q 2 . This is similar to what 
was observed for the pion form factor. The effect of the running of the coupling will be 
discussed further on when we also include the intrinsic transverse momentum dependence. 
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This will allow us to adapt the renormalization scale entering the running coupling to the 



gluon virtuality, as discussed in section [IV B 



By comparing the different vector meson channels in Fig. [16], one sees that the p° L channel 
yields the largest cross section. The u>l channel in the valence region (xb ~ 0.3) is about a 
factor of 5 smaller than the p° L channel, which is to be compared with the ratio at small xb 
(in the diffractive regime ) where p° : uj = 9 : 1. The p\ channel, which is sensitive to the 
isovector combination of the unpolarized OFPD's, yields a cross section comparable to the 
ujl channel. The p\ channel is interesting as there is no competing diffractive contribution, 
and therefore allows to test directly the quark OFPD's. The three vector meson channels 
(p° L , p~£, ujl) are highly complementary in order to perform a flavor separation of the unpo- 
larized OFPD's H u and H d . 

For the pseudoscalar mesons which involve the polarized OFPD's, one again remarks in 
Fig. |l^ the prominent contribution of the charged pion pole to the 7r + cross section. For 
the contribution proportional to the OFPD H, it is also seen that the 7r° channel is about 
a factor of 5 below the ir + channel due to isospin factors. In the ir° channel, the u- and 
(i-quark polarized OFPD's enter with the same sign, whereas in the ir + channel, they enter 
with opposite signs. As the polarized OFPD's are constructed here from the corresponding 
polarized parton distributions, the difference of our predictions for the tt° and tt + channels 
results from the fact that the polarized ci-quark distribution is opposite in sign to the po- 
larized w-quark distribution. For the rj channel, the ansatz for the OFPD H based on the 
polarized quark distributions yields a prediction comparable to the n° cross section. 

For the DVCS, the leading order amplitude is constant in Q and is predominantly trans- 
verse. Therefore, the L.O. DVCS transverse cross section dax/dt shows a 1/Q 4 behavior 



in Fig. |16[ To test this scaling behavior, one needs of course a kinematical situation where 
the DVCS dominates over the BH. As the L.O. DVCS amplitude does involve hard gluon 
exchange as for meson electroproduction, it is a rather clean observable to study the onset of 
the scaling in Q 2 . In Fig. [T7|, we compare the L.O. DVCS transverse cross section (multiplied 
with the scaling factor Q A ) with the result including the transverse momentum dependence 



in the handbag diagrams as described in section IV A. One observes the onset of the scaling 



as Q 2 runs through the range 2-10 GeV 2 , similar to what was also observed for the 7r°7*7 
FF, which is also described at leading order by handbag type diagrams (see Fig. ([7])). It 
will therefore be interesting to measure electroproduction reactions in the range Q 2 ~ 2 - 
10 GeV 2 to study the onset of this scaling. In addition, the preasymptotic effects (in the 
valence region) can teach us about the quark's intrinsic transverse momentum dependence. 

In Fig. [18|, we show how the power corrections modify the L.O. prediction for the longitu- 
dinal p° L electroproduction. One sees that the inclusion of the intrinsic transverse momentum 
dependence leads to an appreciable reduction of the cross section at the lower Q 2 values, 
before the scaling regime is reached. Therefore, the question arises as to how important are 
the other competing mechanisms in this lower /intermediate Q 2 region. We show in Fig. |18 
the estimate of the soft overlap mechanism of Fig. [12|, where the meson is produced without 
invoking a gluon exchange. As discussed in section [IV C| , in this work we are only able to 
estimate the soft overlap contribution from the region \x\ < £ . To estimate the contribution 
which is neglected (|x| > £), we show predictions for different values of xb, as for larger 
values oi xb, the region which is neglected becomes smaller (note on Fig. [18| that at larger 
Xb, the minimal Q 2 value which is kinematically possible in order to be above threshold, 
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increases). One sees from Fig. [18| that the overlap contribution to the longitudinal electro- 
production amplitude drops approximately as 1/Q S , which is indeed the expected result at 
large Q 2 j4|. At i B = 0.3 and Q 2 ~ 3 GeV 2 , our estimate for the soft overlap contribution 
is already more than a decade below the hard electroproduction amplitude including the 
transverse momentum dependence. At xb = 0.6, where the approximation that we make 
in the calculation of the soft overlap contribution is on much safer side, one still observes 
that the soft overlap contribution does not dominate over the hard amplitude and is more 
than a factor of 5 below the hard amplitude including the /c^-dependence, as one approaches 
Q 2 « 10 GeV 2 . This behavior is quite different compared with the overlap contribution to 
the pion electromagnetic FF. From Fig. [11], we indeed see that a similar calculation (using 
also a gaussian ansatz for the fcj_-dependence, with the same parameter a) yields at Q 2 ~ 
10 GeV 2 an overlap contribution to the pion electromagnetic FF which is twice as large as the 
leading order PQCD prediction including the /c^-dependence. Our estimate for the overlap 
contribution indicates therefore that one is in a more favorable situation with longitudinal 
meson electroproduction compared with the pion electromagnetic FF case. Numerically, 
the main reason for this is that in the FF case, the gluon in the L.O. diagram is purely 
spacelike yielding a real quantity whereas in the L.O. meson electroproduction diagram, the 
exchanged gluon can be both spacelike or timelike according to the value of the quark's mo- 
mentum fraction x, yielding a complex amplitude. The form factor is therefore comparable 
to the real part of the hard electroproduction amplitude. As the imaginary part of the hard 
electroproduction amplitude is numerically by far larger than its real part in the kinematics 
considered here, the competing (real) soft overlap mechanism might well dominate the real 
part of the amplitude at the lower Q 2 , without dominating the predictions calculated with 
the total L.O. electroproduction amplitude. 

At present, no experimental data for the p° L electroproduction at larger Q 2 exist in the 
valence region (xb ~ 0.3). However, at smaller values of xb, the reaction 7* p — > p° L p 
has been measured at rather large Q 2 . We will therefore compare our results to see how 
these data approach the valence region, where one is sensitive to the quark OFPD's. For 
the purpose of this discussion, we call the mechanism of Fig. |] which proceeds through the 
quark OFPD's, the Quark Exchange Mechanism (QEM). It is well known that, besides the 
QEM of Fig. [|, p° electroproduction (and, more generally, neutral vector meson electropro- 
duction) can also proceed through a two-gluon exchange mechanism, where it has mostly 



been studied in the past. For this latter, Frankfurt, Koepf and Strikman calculated lon- 
gitudinal p° L electroproduction at low xb and large Q 2 through a Perturbative Two Gluon 
Exchange Mechanism (PTGEM). In this PTGEM, the two quarks connecting to the lower 
blob in Fig. (^) are replaced by two gluon lines. We have implemented this model for the 
PTGEM, in which the longitudinal forward differential p° L electroproduction cross section 
was found as 



iirr L 

' rN ^ VN 1 ■ U^r v ^ +e -M v a 2 s (Q) V 2 v / (a em Q*N 2 ) 



dt 



\t=o 



x |(l+^^)xG T (x,Q 2 )| 2 T(Q 2 ) , (115) 



where Gt(x, Q 2 ) is the unpolarized gluon distribution and 
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_ 1 fdz<Pk ± 9y{z,k ± )/(z(l-z)) 

2 J dzd 2 k± V v (z, k±) 



is equal to 3 if one uses the normalization of Eq. ( |72|) for the vector meson wavefuntion 
*&v(z, k±) with the asymptotic distribution amplitude. The factor T{Q 2 ) in Eq. (|115|) ac- 



counts for preasymptotic effects and was calculated in Ref. [4c] as : 



t(q2) = ( Ipdz J Q d 2 k ± * v (z,k ± )(-\A ± ) A(z,k ± ,Q*, 
V tidztf 2 d*k ± * v (z,k ± )/(z(l-z)) 

where Aj_ is the two-dimensional Laplacian operator in transverse momentum space, and 
where we introduced the notation 

A(z,k ± ,Q 2 ) = ^ 9L . (118) 

Q2 + (kl + m 2 ) I (z(l-z)) 



The correction of Eq. ( |117| ) is due to the Fermi motion of the quarks in the vector meson and 
is equivalent to the intrinsic transverse momentum degree of freedom that we introduced 
in section [TV]. Likewise, it results in a significant suppression of the cross section as Q 2 
decreases {T{Q 2 ) — ► 1 as Q 2 —* oo). 

Intuitively, it is clear that the PTGEM should dominate at high cm. energies, that is 
small xb, where the gluon distribution dominates. The QEM mechanism, which is propor- 
tional to the quark distributions will dominate in the valence region. Due to the sea quark 
distribution, it also contributes and rises in the large W region but quite less than the PT- 
GEM. In Fig. |19, the forward differential longitudinal p° L electroproduction cross section is 
shown as a function of the cm. energy W for three values of Q 2 (5.6, 9 and 27 GeV 2 ). The 
figure shows the contributions of both mechanisms. The PTGEM and QEM cross sections 



are both calculated with the MRST98 parton distribution parametrizations |T7|], including 



their evolution. For the QEM, we include the fcj_-dependence effects as discussed in section 
P"V| . One sees from Fig. |H] that the PTGEM explains well the fast increase at high energy 
of the cross section but it is clear that it substantially underestimates the data at lower en- 
ergies (around W ~ 10 GeV). This is where the QEM is expected to contribute since xb is 
then in the valence region. The QEM describes well the change of behavior of the data (like 
a plateau) at lower W. The incoherent sum of both mechanisms is also indicated. Fig. [20] 
shows this incoherent sum of the two mechanisms for other Q 2 values, up to the largest 
Q 2 value of 27 GeV 2 , measured at ZEUS. It is seen that both W and Q 2 dependences are 
rather well reproduced. This provides a good indication that the deviation from the PT- 
GEM of the data at lower energies can be attributed to the onset of the QEM. The present 
calculation strengthens our previous conclusion about the onset of the QEM at lower W 
and at large Q 2 . The calculation uses now a phenomenological ^-dependent ansatz for the 
OFPD's based on the most recent MRST98 parametrization for the parton distributions. 
Furthermore, the calculation also includes the power corrections due the intrinsic transverse 
momentum dependence both in the QEM as discussed in section [TV], and in the PTGEM 
using the formalism (Eq. ( |115|) ) of Ref. . 



Finally we show in Fig. ^1] how our calculations for 72 p — >• n + n compare with the few 
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data at larger Q 2 . For the sake of comparison with the calculations of Ref. p9| , we also 



show the L.O. predictions for H and for the pion pole using a running coupling constant 
(extrapolated to the lower Q 2 values) with Aq C d = 0.2 GeV. Our results are compatible 
with those of Ref. [29j where it was also found that in the valence region, the pion pole 
dominates the p — ► n + n cross section. We furthermore show in Fig. |21] our predictions 
including the power corrections. For the power corrections to the pion pole contribution we 
use Eq. (B3), where we include in the pion electromagnetic FF the transverse momentum 



dependence and the soft overlap contribution as shown in Fig. |TT| It is seen that in the 
forward kinematics of Fig. |2T|, these power corrections enhance the leading order predictions 
for the cross sections substantially, which is also indicated by the sparse data. It will be 
most valuable to have more accurate data for the longitudinal cross sections for pion elec- 
troproduction to check the consistency of the power correction calculations. Of particular 
interest should be the ratio of 7r + versus it , because the 7r + process, which is dominated 
by the pion pole contribution, has an amplitude with a large real part in contrast to the ir° 
channel, where the pion pole contribution is absent. 



VI. CONCLUSION 

We have given in this work predictions for leading order observables for DVCS and various 
meson electroproduction reactions in the valence region at large Q 2 , using a ^-dependent 
ansatz for the OFPD's. We have indicated some observables and kinematical conditions 
where experiments are already planned or can be performed. As these exclusive experiments 
can currently only be performed at not too high values of Q 2 , we have estimated the power 
corrections due to the intrinsic transverse momentum dependence of the partons in the DVCS 
and meson electroproduction amplitudes. We have taken the 7r°7*7 transition FF and the 
pion electromagnetic FF as our guidance to estimate these power corrections to both DVCS 
and hard meson electroproduction reactions respectively. In this way, we have generalized the 
skewed parton distribution formalism to include the parton intrinsic transverse momentum 
dependence. For the meson electroproduction amplitude, we have estimated, in addition, 
the competing soft overlap mechanism, which - in contrast to the leading order perturbative 
mechanism - does not proceed through one-gluon exchange. The soft overlap contribution is 
found not to be dominant in contrast to the pion FF case. Our estimates for these different 
power corrections show that by measuring exclusive electroproduction reactions in the range 
Q 2 m 2 - 10 GeV 2 , one will be able to study the onset of the predicted scaling. In addition, 
the preasymptotic effects (in the valence region) can teach us about the quark's intrinsic 
transverse momentum dependence in the nucleon. 

In conclusion, we believe that a broad new physics program, i.e. the study of exclusive 
reactions at large Q 2 in the valence region, where the quark exchange mechanism dominates, 
opens up. Although these exclusive experiments at large Q 2 are quite demanding, we think 
that both the scaling region and the onset of the scaling region promise to be sufficiently 
rich to motivate an extensive experimental investigation at different facilities. 
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TABLES 





dxx H 


f\dxx(H + E) 


SMC 




u v 


.28 


.51 


An, = .77 ± .10 ± .08 


-.13 


(ly 


.11 


-.12 


Ad v = -.52 ± .14 ± .09 


.20 


Uy ~\~ dy 


.39 


.51 


Au v + Ad v =.25 


.07 


U + U 


.34 


.62 


Au v + 2A-u=.79 


-.09 


d + d 


.18 


-.19 


Ad v + 2Ad=-.50 


.16 


{u + u)+{d + d) 


.52 


.43 


E=-29 


.07 



TABLE I. The evaluation of the spin sum rule using the ^-independent factorized ansatz for 
the OFPD's as described in the text. The column denoted by SMC contains the experimentally 
determined [24| contributions of the different quark flavors to the nucleon spin, measured through 
semi-inclusive spin asymmetries. Remark that the model calculation implies a fraction of 0.57 for 
the gluon contribution to the nucleon spin. 
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(a) (b) 
FIG. 2. "Handbag" diagrams for DVCS. 




p p' = p + A 

FIG. 3. Diagram to denote the arguments x, y and A 2 of the double distributions. 
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FIG. 5. f dependence of the OFPD H d at t = using the ansatz (based on the MRST98 
quark distributions) as described in the text. Upper panel : valence down quark OFPD, lower 
panel : total down quark OFPD. The thin lines (£ = 0) correspond with the ordinary (i-quark 
distributions (MRST98 parametrization at Q 2 = 2 GeV 2 ). 
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FIG. 6. Leading order diagrams to hard meson electroproduction. 
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(q 2 = -Q 2 < 0) 




q' (z, -k ± ) 



(q 12 = 0) 

(°) (6) 
FIG. 7. Leading order direct (a) and crossed (b) diagrams to the 7r°7*7 transition form factor. 
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FIG. 8. Results for the 7r°7*7 transition form factor of the leading order PQCD prediction 
(dashed-dotted line) compared with the prediction including transverse momentum dependence 
(full line). Data are from CELLO || and CLEO 
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FIG. 9. Impact parameter (b) dependence of the OFPD H v for the valence down quark, usinj; 
the model ansatz (based on the MRST98 quark distributions) as described in the text. 
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FIG. 10. Leading order direct (a) and crossed (b) diagrams to the tt electromagnetic form 
factor involving one-gluon exchange. Soft overlap contribution (c) to the tt form factor. There are 
analogous diagrams where the virtual photon couples to the lower quark line, which are not shown 
here. 
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FIG. 11. Results for the ir electromagnetic form factor of the leading order PQCD prediction 
without (dotted line) and with (dashed-dotted line ) inclusion of the corrections due to intrinsic 
transverse momentum dependence. The dashed curve shows the result for the soft overlap contri- 
bution of Fig. |lO| c and the total result (full line) is the sum of the dashed and dashed-dotted lines. 
The "data" points from Refs. H4|] (CEA) and E3] (Cornell) were obtained indirectly from pion 
electroproduction. 
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FIG. 13. Comparison between the angular dependence of the leading order predictions for the 
p\ (full lines), ir + (thick, upper dashed lines) leptoproduction and DVCS (thick dotted lines) 
in-plane cross sections at Q 2 = 2.5 GeV 2 , xb = 0.3 and for different beam energies : = 200 
GeV (COMPASS), E e+ = 27 GeV (HERMES), E e - = 6 GeV (JLab). The BH (thin dotted lines) 
and total 7 (dashed-dotted lines) cross sections are also shown. For the 7r + , the result excluding 
the pion pole is also shown (thin, lower dashed lines). 
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FIG. 14. Angular dependence of the leading order amplitudes for p® L leptoproduction at E^ 
= 200 GeV, Q 2 = 2.5 GeV 2 , and for different values of xg. The prediction with a ^-dependent 
ansatz for the OFPD's (full lines) are compared with a ^-independent ansatz for the OFPD's 
(dashed-dotted lines). 



52 



e" + p -> e" + p + y (O = 120°) 




© y (deg) 

FIG. 15. DVCS at JLab for out-of-plane kinematics (4> = 120°). Upper part shows the dif- 
ferential cross section : DVCS including the 7r° pole contribution (dashed line), BH (dotted line), 
total 7 excluding the ir° pole contribution (dashed-dotted line) and total 7 including the 7r° pole 
contribution (full line). In the lower part, the electron single spin asymmetry is shown for the L.O. 
DVCS amplitude (calculated in radiative gauge) excluding the 7r° pole (dashed line) and including 
the 7T° pole (full line). We also show the result for the gauge invariant DVCS amplitude, excluding 
the 7T° pole (dashed-dotted line). 
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FIG. 16. Scaling behavior of the leading order predictions for the forward differential electro- 
production cross section dai/dt on the proton, for vector mesons (left panel) and pseudoscalar 
mesons (right panel) using a coupling constant frozen at a scale 1 GeV 2 . For the 7r + channel, the 
pion pole contribution (full line, (7r + ) po ; e ) is shown separately from the H contribution (dashed 
line, (vr + )„. p .). Also shown is the scaling behavior of the forward transverse cross section dar/dt 
for the leading order DVCS cross section (dashed-dotted line in left panel). 
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FIG. 17. DVCS transverse cross section dar/dt (multiplied with the scaling factor Q 4 ) in the 
forward direction (t = t m i n ). The leading order result is given by the dashed-dotted line and the 
result including the intrinsic transverse momentum dependence is given by the full line. 
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FIG. 18. Longitudinal differential cross section dai/dt in the forward direction (t = t m in) for 
for ^)* L p — ► p^p . The leading order result is given by the dashed line and the result including 
the intrinsic transverse momentum dependence is given by the full line. The calculation of the soft 
overlap contribution is given by the dashed-dotted line. 
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FIG. 19. Longitudinal forward differential cross section for p a L electroproduction. Calculations 
compare the quark exchange mechanism (dotted lines) with the two-gluon exchange mechanism 
(dashed lines) and the sum of both (full lines). Both calculations include the corrections due to 
intrinsic transverse momentum dependence as described in the text. The data are from NMC 
(triangles) p9[ , E665 (solid circles) [ pO[j , ZEUS 93 (open circles) |5l| and ZEUS 95 (open squares) 
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y* + p -> p° + p : Q 2 = 5.6, 9, 13, 17, 27 GeV 2 




FIG. 20. Longitudinal forward differential cross section for p a L electroproduction. Calculations 
show the sum of the quark exchange and two-gluon exchange mechanisms for different values of 
Q 2 (in GeV 2 ) as indicated on the curves. The data are from NMC @, E665 fU, ZEUS 93 || 
and ZEUS 95 
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FIG. 21. Q 2 -dependence for the longitudinal forward differential electroproduction cross section 
daL/dt for — > 7r + n at xb = 0.24. The leading order pion pole contribution (dotted line) and H 
contribution (dashed-dotted line) are shown using a running coupling constant with Aqcd = 0.2 
GeV. Also shown (dashed line) is the calculation including the power corrections to the pion pole 
contribution (intrinsic transverse momentum dependence and overlap contribution to the pion FF), 
as well as the sum of the pion pole and H contribution including the power corrections (full line). 
The data are from Ref. p5l : the points at Q 2 = 1.2 and 2.0 GeV 2 correspond with xb ~ 0.24, 
whereas the point at Q 2 = 3.3 GeV 2 corresponds with xb ~ 0.35. 
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